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£S . ABSTRACT 

The status of the calculation of next-to-leading-order QCD corrections to 
hadronic final states in deeply inelastic scattering is reviewed, and an overview of 
the phenomenology (including the measurement of the strong coupling constant 
I . and the gluon density via jet rates) is given. We also describe a new universal 

J> ' program (DISASTER++) for the calculation of (2+l)-jet observables. 

(N : 

^ ■ 1. Introduction 

oo : 

The electron-proton collider HERA has started to accumulate large luminosity. 
Q\ ■ This opens up the possibility to study the hadronic final state in deeply inelastic 

scattering to high precision. The motive for this is twofold: QCD can be tested at large 
scales in the spacelike regime, and an extraction of physical quantities, in particular 
Qh! the running strong coupling constant a s (Q 2 ) and the gluon density f g ((,,Q 2 ), is 

possible. The aim of this paper is to give an overview of the status of next-to-leading- 
order calculations both for jet-like quantities and for one-particle-inclusive processes. 
We will also discuss the phenomenology of processes with (2+1) jets a in the final 
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$_i ■ state, and of inclusive particle spectra. 

In the last few years, the technology for the calculation of QCD corrections in 
next-to-leading order of perturbation theory has developed considerably. There are 
explicit algorithms available which permit the calculation to be done in a "universal" 
way: the infrared singularities are subtracted in such a way that arbitrary infrared- 
safe observables can be calculated numerically. In principle, all existing algorithms 
are variations on a common theme, namely the interplay of the factorization theorems 
of perturbative QCD and the infrared-safety of the observables under consideration. 
This will be reviewed in Section 0. 

The practical implementations of this principle can be quite different. In Section |3|, 
we will discuss the two possible ways to do this: the phase-space-slicing method and 
the subtraction method. The extraction of the singular terms can be done in different 
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a Here "(n + 1)" stands for n hard jets and the proton remnant jet. 



ways, and we will describe in some detail a new calculation based on the subtraction 
formalism and on a general partial fractions formula. 

An overview of applications for (2+l)-jet-like observables will be given in Sec- 
tion [|, where we state some results for the cut dependence of jet cross sections and 
for transverse momentum spectra. One-particle-inclusive processes will be discussed 
in Section || The theoretical predictions in the framework of the fragmentation func- 
tion picture fail to describe the experimental data for small values of the photon 
virtuality Q 2 and for small values of the energy of the observed hadron, although 
there is excellent agreement at large Q 2 and large E^. We will give a semi-quantitative 
estimate of the region where the fragmentation function picture should be applicable. 
The paper closes with a summary and an outlook. 

2. Infrared-safe Observables and Factorization in QCD 

Perturbative QCD permits the calculation of processes with partons in the final 
state. The gap to experimental data with final-state hadrons is bridged by means of 
infrared-safe observables O. Examples are jet cross sections for various jet definitions, 
and event shape variables such as Thrust. The expectation value of an observable O 
for experimental data is given by 

i N 
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where N is the number of events, nj is the number of final-state hadrons of event /, 
and the hf^ are the momentum vectors of the hadrons. This is to be compared with 
the expectation value for parton final states: 

(0) th = W dPSW <r<»> (pi, . . . ,p n ) W (Pi, • • • ,Pn) ■ (2) 

Here n is the number of final-state partons with momenta pi, p n , dPS^ is 
the phase space measure, and a<ri is the hard scattering cross section for partons 
in the final state, calculated in perturbation theory 6 . For the next-to- leading-order 
calculation of the (2+l)-jet-production cross section to be described later on, the 
relevant values are n = 2 for the Born term and the virtual corrections, and n = 3 for 
the real corrections. Both virtual and real corrections exhibit infrared singularities. 
The soft singularities and the collinear singularities in the final state cancel, owing to 
the Kinoshita-Lee-Nauenberg mechanism. The collinear singularities from the initial 
state can be consistently absorbed into redefined parton densities. The resulting hard 
scattering cross section is finite. 

6 We have not indicated technical complications such as different types of final-state partons (quarks, 
gluons) and the convolution with parton densities if there are hadrons in the initial state. 



The factorization theorems of perturbative QCD (see Ref. |]J and references 
therein) show that parton cross sections a have a very simple behaviour in collinear 
and soft phase space regions: 



In the collinear limit, for the partons with labels j and k being collinear, the 
limit for o is 

O ~ 77^ Pkj^iiu) CTBorn, (3) 

Ztt s jk 

where u is the momentum fraction of the parent parton % carried by parton j, 
Sjk = 2pjPk and (XBorn is the Born cross section corresponding to the real cor- 
rection a in the limit under consideration. The Pkj^i{u) are the unsubtracted 
Altarelli-Parisi splitting functions (see, for example, Ref. ||). There is a slight 
technical complication: factorization in this form holds only for polarized cross 
sections. The sum over the parton polarizations introduces residual azimuthal 
correlations, even in the collinear limit. These have to be taken into account in 
the construction of the subtraction terms, see below. After integration over the 
azimuthal angle, factorization of the form of Eq. ([3D holds also in the unpolar- 
ized case. The splitting functions Pkj^iiu) are universal, process-independent 
functions. This has the important consequence that physical parton densities 
and fragmentation functions can be defined in a process-independent way. 

In the soft limit, the situation is slightly more complicated. For a gluon c with 
label k becoming soft, the factorization formula now reads 



The sum runs over all pairs of partons {i, j}, and the are constant coeffi- 
cients. The structure of this formula can be easily understood in terms of an 
eikonal approximation, where the matrix element Ai factorizes as 

M-^AW (5) 

PiPk 

The sum over the gluon polarization A, e(A) being the gluon polarization vector, 
leads to the form given in Eq. (|4]). 

In order not to spoil the cancellation of soft and collinear singularities from the real 
corrections against those of the virtual corrections, the observables O have to fulfill 
certain conditions. This is the topic we now turn to. 



c Soft quarks do not lead to a soft singularity. The matrix element is singular as 1/E, but this 
singularity is compensated by a factor E in the phase space volume. 



An observable O is called infrared-safe, if the functions O^ (pi, . . . ,p n ) of the 
parton momenta respect soft and collinear limits: 

0^(p u ..., Pi ,...,p n ) — O^Cpx,...,^,...,^), (6) 

P<-*0 

C (n) (pi, . . . . . . . .,p n ) — -> cX™" 1 ) (p Xj . . . 7 p h . . . . . . ,p n ,pi +pj) . 

Pi\\pj 

Momenta denoted by p are to be omitted. The property of infrared safety for observ- 
ables has the consequence that the factorization from Eqs. (|3]) and (|4]) also works if a 
is replaced by a O. This ensures the cancellation of the infrared singularities for the 
case of a convolution of the parton cross section with an observable, as in Eq. @. 

3. Universal Calculations in Next-to-Leading Order for Jet Quantities 

The factorization theorems of perturbative QCD make it possible to perform cal- 
culations for a particular process in such a way that arbitrary infrared-safe observables 
can be evaluated numerically. For this it is necessary to extract and cancel the infrared 
singularities in an observable-independent way. Technically, there are two different 
procedures to achieve this: the phase-space-slicing method and the subtraction 
method These techniques can be illustrated by means of a simple example 
Assume that the integral 

/ = / dxx~ e -f{x) (7) 
Jo X 

is to be evaluated. The x-integration stands for the phase-space integral, the factor 
x~ e is the regulator 6 , the term 1/x is the infrared singularity (a "propagator term"), 
and the integrable function f(x) contains regular parts of the cross section and of the 
observable. 

• In the case of the phase-space-slicing method, the range of integration is split 
into two parts [0, a] and [a, A] by means of an arbitrary small technical cut-off 
parameter a. For the lower interval, the function f(x) can be approximated by 
/(0) up to terms 0(x), which regulate the 1/x singularity, and, after integration, 
give rise to terms 0(a). For the upper interval, the regulator e can be set to 
zero. The final result for the integral is thus 

I = --/(0) + /(0) lna + / dx - f(x) + 0(a) + 0(e). (8) 

6 J a X 

The first term is the singular part, which would cancel against a similar sin- 
gular part of opposite sign from the virtual corrections. The remaining terms 

d Reference || actually contains the first fully general algorithm to calculate subtraction terms. 
e The calculations are performed in d = 4 — 2e space-time dimensions, so that UV and IR singularities 
are regulated and made explicit in the form of poles in e. Infrared singularities are regulated for 
e < 0. 



are finite. The integral in the third term can be evaluated numerically. The 
method is exact only in the limit a — * (in practice, very small values for a are 
sufficient). The large logarithm in the second term is compensated by a similar 
logarithm from the numerical integration in the third term. 

The subtraction method requires the definition of a subtraction term, such that 
a splitting of the integral into a finite and a singular part is possible. As an 
example we give a specific way to define the subtraction term: 

A dxx~ e - (f(x) - f(0)) + [ A dxx~ e - /(0). (9) 



10 X Jo X 

The subtraction ("local in phase space") in the first integral renders the inte- 
grand integrable for e — > 0, and the final result is thus 

J=--/(0) + /(0)lnA + [ A dx- (/( Z )-/(0))+O(e). (10) 
e Jo x 

The integral in the third term can again be performed numerically. Because 
/(0), the subtracted term in the integrand, is the residue of f(x)/x at x — 0, 
this method is also called the residue method. 

The advantage of the phase-space-slicing method is that it can be implemented in 
a straightforward way: all that is required is the integral of the real corrections 
over small cones in phase space (corresponding to the approximation x ~ 0) and a 
numerical implementation in the form of a Monte-Carlo program of the cross section 
itself. The disadvantage is that the cancellation of the term proportional to In a is 
delicate; the statistical fluctuations require a large number of Monte-Carlo events. 
Moreover, the convergence of the final result in the limit a — > has to be checked, in 
principle, for every calculated expectation value. These technical problems are not 
present in the case of the subtraction method, because no small cut-off has to be 
used-^. However, the implementation of the method requires the construction of the 
subtraction term, which is a non-trivial task. If this can be afforded, the subtraction 
method is the method of choice. 

For the particular case of (1+1)- and (2+l)-jet production in deeply inelastic 
scattering, there are by now several calculations in the form of weighted Monte-Carlo 
programs available: 

• PRO JET |J: The jet definition is restricted to the modified JADE scheme; the 
program is based on the calculation in Refs. 



DISJET [p70 1 : Again the jet definition is restricted to the modified JADE scheme; 



the program is based on the calculation in Refs. [11, 12 . 



^Strictly speaking, a tiny cut-off of the order of to 1CT 10 has to be introduced in order to make 
the cancellation (f(x) — /(0)) fx work, because there is only a finite number of significant digits 
available. 



MEPJET ||13|| : This is a program for the calculation of arbitrary observables which 
uses the phase-space-slicing method. The calculation |n| uses the Giele-Glover 
formalism Jl5| for the analytical calculation of the IR-singular integrals of the 



real corrections, and the crossing-function technique [16j to handle initial-state 
singularities. The latter requires the calculation of "crossing functions" for each 
set of parton densities. 

D I SENT : This 

program is based on the subtraction method. The subtraction 



term is defined by means of the dipole formalism 9 JL8], |T 



• DISASTER++ ||20|| : This is a C++ class library' 1 . The subtraction method is em- 
ployed, and the construction of the subtraction term resembles the method of 
Ref. ||, i.e. it is obtained by the evaluation of the residues of the cross section in 
the soft and collinear limits. Double counting of soft and collinear singularities 
is avoided by means of a fully general partial fractions method. 

Why a new calculation? There are two reasons: (a) The existing programs have 
the restriction that the number of flavours is fixed (Nf = 5 in the case of MEPJET and 
Nf fixed, but arbitrary for DISENT). For studies of the scale-dependence it is necessary 
to have a variable number of flavours, in order to be consistent with the scale evolution 
of the strong coupling constant and the parton densities. DISASTER++ makes the Nf 
dependence explicit in the "user routine" on an event-by-event basis, and thus results 
for arbitrary renormalization and factorization scales can be binned simultaneously, 
(b) DISASTER++ is already set up such that the extension to one-particle-inclusive 
processes will be possible without the necessity of re-coding the contributions which 
are already present for the jet-type observables. 

In the following we will illustrate how the subtraction term for the real corrections 
is constructed. We will be concerned with infrared-safe observables that are zero if 
more than one parton is soft, or if more than two partons are collinear, or combinations 
thereof. These are the most relevant ones in practical applications. In principle, the 
procedure is simple. In energy and angle variables, soft singularities are of the form 
l/Ef, for parton i being soft, and collinear singularities are given by terms 
for partons i, j being collinear, where = (1 — cos$ij)/2, and i?^- is the angle 
between the two partons. The singularities can be extracted by performing the limits 
Ei — > and — > of the terms Efa, Vij<T and Efvija, and the subtraction terms 



9 The subtraction term is written as a sum over dipoles (an "emitter" formed from two of the original 
partons and a "spectator" parton). Besides the factorization theorems of perturbative QCD, the 
main ingredient is an exact factorization formula for the three-particle phase space, which allows for 
a smooth mapping of an arbitrary 3-parton configuration onto the various singular contributions. 
h The acronym stands for "Deeply Inelastic Scattering: All Subtractions Through Evaluated 
Residues". Most of the program is written in C++. A FORTRAN interface is in preparation; thus 
there will not be any problem to interface the class library to existing FORTRAN code. 



will consequently be given by 
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Ef a 



E i v i: j a 



There is, however, the problem of overlapping singularities, or double counting. Con- 
sider a term of the form \j{E\ v\2 v 13). Here the problem is that if parton 1 is soft, then 
there are two regions which give rise to an additional collinear singularity: v 12 = 
and v 13 = 0. In order to include every possible singular configuration only once in the 
subtraction term, it is convenient to separate the singular terms. This can be done 
in two ways: (a) division of the phase space, or (j3) separation of the singularities 
in the matrix element. In Ref. || a combination is used: 1/(^12^13) is written in 
terms of partial fractions as l/[fi2(fi2 + v n)] + 1/^13(^12 + ^13)]) an d the energy of 
the i^-integration is restricted to be smaller than E 2 /2. In principle, phase-space 
cuts of this kind do not pose a problem in Monte-Carlo programs. However, it is 
preferable that the integrand be a smooth function. We achieve this by using the 
following general formula for partial fractions: 



X\ X2 ' ' ' X^i Q-E.S (^""l *^'(T2 ) ' ' ' (*^<Tl ~t~ • • • ~t~ ^/J n 



(12) 



The sum runs over all n! permutations of n objects 1 . The most straightforward 
way to apply this formula to the case at hand [(2+l)-jet production] would be to 
set the set of variables . . . , x 9 } to {Ef, E%, E\, v i, v 02 , v 03 , Vu, v 13, t^} (po is the 
momentum of the incident parton), where the energies are rescaled such that they 
are dimensionless. The product of the cross section and the observable a O is then 
rewritten as [xi ■ ■ • x n a 0\j{x\ ■ ■ -x n ), and the partial fractions identity (|T^) is ap- 
plied to the denominator. The numerator is regular: the two-parton singularities are 
regulated by the x iy and the remaining singularities are regulated by the observable 
(recall that we restrict the discussion to this kind of observables) . The general par- 
tial fractions formula now gives us a "hierarchy" of singularities: the "leading" one 
l/x ai , the "subleading" one l/(x CTl + x a2 ), and so forth. It turns out that for the 
observables under consideration, we have to consider only the leading and subleading 
singularities. Let us write this singular part, for a specific term in the sum, as k^, and 
the remaining terms from the product oO as ta', the index A parametrizes the set 
of permutations. We then have a O = J2a^a^a- Each of the terms is singular only 
if the "leading variable" is zero. The subtraction term can therefore be constructed 
with respect to the leading and subleading variable in fc^. More precisely this is done 
in the following way. For each index A there is a particular parton label ia related 

' This decomposition has a nice technical property. Holding o\ fixed, excluding l/x ai from the sum, 
and setting x ai to zero in the remainder, the sum over the restricted set of permutations yields 
1/ {x\ ■ ■ ■ x ai ■ ■ ■ x n ). This property is useful to recombine terms after the soft and collinear limits 
have been performed. 



to it (for example, for l/[Ef (Ef + Di 2 )] or l/[vi 2 (Ef + v\ 2 )] this would be %a = 1)- 
The phase space dPS (n) is factorized according to dPS iA dPS (n - 1} , i.e. the one-parton 
phase space corresponding to parton %a is pulled out of the full phase space. Finally 
the phase space integration from Eqn. (El) is rewritten as 



J dPS (n) a O = E J dPS *A k A (j dPS (n_1) rA - J dPS^ 1 
+ E/ dPS iA fc A [| dPS^- 1 )^ 



TA 



soft/coll, limit j 



soft/coll, limit 



(13) 



The first integral is finite and can be calculated numerically. The second integral 
contains all infrared singularities. The term in the square bracket is simple because 
of the factorization theorems of QCD, and the one-particle integral over the kernel 
Ha and the factorization contribution from the term in the square brackets can be 
performed easily- 7 . 



4. Numerical Results and Applications: ct s and the Gluon Density 

We now turn to some applications of (2+l)-jet observables. Loosely speaking, 
observed jets of hadrons can be identified with partons at hard scales. The sub- 
sequent fragmentation process is assumed to be sufficiently soft in order to keep 
non-perturbative effects small. Jets cannot be defined in a canonical way, rather they 
are objects by definition. It is convenient to define jets in terms of iterative cluster 
algorithms^. There are three ingredients in the definition of a cluster algorithm: (i) a 
distant measure which determines the relative distance of two clusters in momentum 
space (for example their invariant mass = (pj + Pj) 2 , (ii) a mass scale M 2 that 
determines whether two clusters are to be combined into a single one (for example if 
Sij < M 2 ), and (iii) a recombination procedure that prescribes how two clusters 
and pj are to be merged into a single cluster p* (for instance p* = Pi+Pj)- There are 



various choices for (i), (ii) and (iii). The example given is the JADE algorithm [21 
in the so-called E-scheme. This particular algorithm is known for large hadronization 
corrections. This problem can be reduced by explicitly forcing clusters to behave as if 
they were massless (JADE and P-schemes). The large corrections partly come from 
the fact that for JADE- type algorithms, soft partons are combined first, even if they 



3 There are a few complications, though. The first one is that the required sum of 9! terms, which has 
to be done numerically, contains too many terms to be efficient (or even possible) . The solution is to 
perform the partial fractions decomposition separately for the energy and angle terms. The second 
complication is the presence of correlations in azimuthal angle for collinear singularities, which have 
been mentioned above. They can be dealt with by introducing a fictitious azimuthal variable for the 
collinear configurations. 

k "Non-cluster" algorithms of the cone-type suffer from the problem that the outcome may depend 
on the choice of a seed in the jet-finding process. This effect is particularly large at small transverse 
momentum. Cluster algorithms are well-defined without any ambiguity in the procedure. 



differ considerably in their direction. Intuitively, such partons should possibly be 
combined with other partons nearby, even if they are hard. This can be achieved by 
means of a distance measure based on relative transverse momentum, leading to the 
kx algorithm, defined in the Breit frame of reference. For deeply inelastic scattering, 
it has been introduced in Ref. W% 1 ■ 



In the case of hadrons in the initial state, the corresponding remnants need a 
special treatment, because partons in the very forward direction lead to initial-state 
singularities. This problem is taken care of by either including the remnant jet in the 
clustering procedure (as in the case of the modified JADE algorithm) or by discarding 
partons which are too close to the remnant (as for the kx algorithm). In any case a 
well-defined prediction for the (2+l)-jet cross section is obtained. 
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Figure 1: Cut- dependence of R^+i) for the modified JADE cluster algorithm (a) and 
for the kx algorithm (b). The jet resolution mass is M 2 = cW 2 for the JADE algo- 
rithm and M 2 = cQ 2 for the kx algorithm. Leading order [ /, next-to-leading 

order: E-scheme [■■■], P-scheme [ ], JADE-scheme [-■-]■ The recombination 

scheme dependence of the kx algorithm is small, therefore only the E-scheme is shown 
in this case. 

In order to reduce systematic errors, it is convenient to consider the (2+l)-jet rate 
R(2+i) = o~2+i/cr tot . Figure |I] shows the cut- dependence of R(2+i) for two different jet 



These algorithms have the property of factorization, which means that the (n+ l)-j et cross section 
factorizes in exactly the same way into a product of a coefficient function and a parton density as 
a structure function does, without any residual explicit x^-dependence. This kind of factorization 
is technically convenient for the calculation of resummed cross sections. We wish to stress that any 
acceptance cut on the final-state jets destroys this factorization property. What is important to 
have in a well-defined perturbative prediction is that the singular parts factorize in a universal way, 
which permits the definition of process-independent parton densities. This is, of course, the case for 
all infrared-safe observables, including the non-factorizing kind of jet algorithms. 



algorithms"? The recombination scheme dependence of the modified JADE algorithm 
is large, as shown in Fig. |I]a. For the E-scheme the QCD corrections may be as large 
as 100 %; they are considerably smaller for the "massless" recombination schemes 71 . 
The perturbative stability of the kr algorithm is better, here the QCD corrections 
are of the order of a few per cent (Fig. [TJb). 
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Figure 2: Transverse momentum distribution in the hadronic CMS for the JADE 
algorithm (a), and in the Breit system for the hp algorithm (b). Thrust distribution 

in the Breit frame (c). Shown are the leading order [ / and the next-to-leading 

order [ J. 



Deeply inelastic scattering is a multi-scale problem. Besides the photon virtuality 
Q 2 there are other hard scales related to the scattering process, for example the 
transverse momentum pt of partons or jets. These hard scales can be very different 
from each other. The transverse momentum distribution of the jets for fixed lepton 
variables of x B = 0.01 and Q 2 = (20 GeV) 2 is shown in Fig. g for the JADE P- 
scheme with a recombination scale of 0.02 W 2 in the hadronic CMS (a) and for the kx 
algorithm with a recombination scale of 0.5 Q 2 in the Breit frame (b). The transverse 
momentum distribution is quite broad, and therefore pj, is not necessarily of the order 

m The centre-of-mass energy is 300 GeV, the photon virtuality is restricted by (10 GeV) 2 < Q 2 < 
(20 GeV) 2 , the total hadronic energy W is assumed to be larger than 70 GeV and the lepton variable y 
is restricted to values smaller than 0.7. The parton densities are MRS A' pj| , and the renormalization 
and factorization scales are set to Q 2 . 

" The large recombination scheme dependence of the JADE algorithm sometimes leads to the con- 
clusion that the algorithm is ill-defined. This is a misconception. Different recombination schemes 
define different observables, and there is no reason why different observables should not lead to 
different predictions. However, the strong dependence of the jet rate on the recombination scheme 
points to the fact that jet masses are actually important. 



of Q 2 . This introduces a theoretical uncertainty, because it is not clear whether Q 2 or 
p\ is the right scale to be used in the renormalized coupling a s {ii 2 ) and in the parton 
densities /i(£>A*/)- As long as no resummed calculation is available, this uncertainty 
has to be considered as a contribution to the systematic error in case that physical 
quantities are extracted by using fixed-order next-to-leading-order calculations . 

As an example for an event shape variable we show a Thrust distribution in Fig. 0c. 
The phase space region is the same as the one of the jet rates shown before. The 
Current Thrust variable is defined by T z = 2Y^iPi\\/Q [HE' where the sum is over all 
partons in the current hemisphere (the hemisphere containing the virtual photon). 
The quantity pq is the longitudinal component of the parton momentum p, along the 
virtual photon direction in the Breit frame of reference. The Thrust variable goes to 
one for a (l+l)-jet-like configuration. Values of T z < 1 are therefore obtained from 
processes of 0{a s ). The QCD corrections are moderate. It turns out that the NLO 



calculations do not describe the experimental data well (see, for example, Ref. p5|). 
The inclusion of a power suppressed term of the form A/Q with a large coefficient 
A brings data and the theoretical prediction in good agreement. A quantitative 
derivation for the coefficient A in terms of an effective coupling constant at low scales 
is given in Ref. |26| . 



Applications of NLO calculations in deeply inelastic scattering include the mea- 
surement of the strong coupling constant |27],|28| and a direct fit of the gluon density 
|[29|| . The a s measurements show a clear evidence for the running of the coupling 
constant as a function of the scale Q 2 . To complement the measurement via jet rates, 
measurements of a s based on event shape variables are currently under way [f23p. 
The gluon density has been determined via jet rates by means of a Mellin transform 
method P0,|5T which makes the repeated evaluation of the NLO cross section for the 



purpose of the fitting procedure feasible 9 . The quark distributions are well known, 
and therefore the quark-initiated contribution can be subtracted. The obtained direct 



°Not shown here because of lack of space is a comparison of the scale dependence in leading and 
next-to-leading order. For the jet cross sections shown in Fig. [I] the scale dependence is reduced by 
about a factor of 2, if the renormalization and factorization scales are varied between 0.5 Q and 2 Q. 
p For the Thrust variable, the requirement of a power-suppressed term A/Q makes an a s measure- 
ment conceptually difficult, because the coefficient A is parametrized by a universal effective coupling 
constant a s e g , whereas the perturbative contribution has the renormalized coupling constant as a 
coefficient. There are matching procedures to disentangle the two terms, but the size (about the 
same order of magnitude as the perturbative contribution) of the non-perturbative contribution is 
disturbing. 

q This method is actually quite general and permits the calculation of arbitrary observables, the only 
limitation being that the factorization scale has to be fixed. 



fit is in good agreement with gluon density parametrizations from global fits r . 



5. One-Particle-Inclusive Processes 



One-particle-inclusive processes are another very promising field. The correspond- 
ing cross sections can be considered to be related to a special type of observable O, 
the fragmentation functions Di (£,//£>). In the terminology used above, they are not 
infrared-safe, so that collinear singularities remain, which do not cancel. However, 
the singularities are of a universal form, and can be absorbed into the fragmenta- 
tion functions. The redefined fragmentation functions are finite and universal, i.e. 
process- independent. 
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Figure 3: Theoretical 'prediction vs. experimental data for the Q 2 -dependence of 
c h /o- tot for various bins of x p . The kinks in the theoretical curves come from the 
overlay of several bins in xb for the same bin in Q 2 . The region within the triangle 
cannot be described by means of the fragmentation function picture; see text. 



Figure |3] shows a comparison of recent preliminary experimental data for charged 
hadron production from the ZEUS collaboration with a next-to-leading order calcu- 
lation. The employed variable is x p = 2Ph/Q, where Ph is the momentum of the 

r For a specific value of the momentum fraction, the error bands show a cross-over. This phenomenon 
is well known from global fits p2| and is a consequence of the limited number of parameters in the 
ansatz for the gluon density. 



observed charged hadron h in the current hemisphere of the Breit frame. The cross 
section o h is a convolution of parton densities, fragmentation functions and pertur- 
batively calculated coefficient functions. The parton densities are from the MRS A' 
set and the fragmentation functions are from Refs. p3|,p4[. The coefficient functions 



have been calculated in Ref. [J35|J. The numerical results have been obtained by us- 



ing the implementation in the program CYCLOPS of a recent recalculation |3q| . The 



comparison of the theoretical prediction with experimental data shows a very good 
agreement, except for small values of Q 2 and small values of x p . 

What is the reason for this discrepancy? It is instructive to consider the situation 
in terms of the rapidity variable yh of the observed hadron. We consider the leading- 
order process, with a single parton fragmenting into the hadronic final state. Under 
the assumption that the hadron has a typical transverse mass my of 0(500 MeV), the 
relation between x p and can be calculated. In the Breit frame, a positive value of yh 
stands for production in the current direction, and a negative value for production 
in the target direction. In order to make the fragmentation function picture a valid 
description, the hadron should be produced closely in rapidity to the parent parton, 
which means that yh should be larger than about 7/ min = 1 units of rapidity. This can 
be translated into a lower bound on x r . 



Q J 1 — (tanh y 



x ■ =— (14) 



mm I 



This excludes the region roughly indicated by the triangle in the figure. The theo- 
retical prediction thus fails in this region because the fragmentation function picture 
is not applicable: the hadron is not produced sufficiently close to the parent parton. 
Another reason is that mass effects (which are not included in the fragmentation 
function formalism) become important at small Q 2 and small x p for x p ~ 2m n /Q. 
This excludes a similar region in the plot. Outside of this critical region, in particular 
at large values of Q 2 , the good agreement shows the universality of fragmentation 
functions. 

A future application might be the determination of a s via scaling violations of the 
fragmentation functions. This has been studied in some detail in Ref. [p7H . It turns 



out that the main uncertainty is the dependence on a parametrization of the parton 
densities. This dependence is reduced at large values of Q 2 , but then the available 
luminosity is the limiting factor. Another possibility might be the measurement of a s 
from the hard scattering matrix element by means of px spectra of charged particles. 
For this the NLO calculation is not yet available; however, as indicated above, the 
implementation of the matrix elements in DISASTER++ is already such that it may 
become available in the near future. The advantage of such a measurement would be 
that the fragmentation process is no longer modelled by means of event generators as 
in the case of jet cross sections or event shape variables, but described by a universal 



parametrization of fragmentation functions with a QCD-predicted scale evolution. 



Summary and Outlook 

We have reviewed the status and the applications of next-to-leading-order calcu- 
lations in deeply inelastic lepton-nucleon scattering. In the last few years, programs 
have become available that permit the calculation of arbitrary infrared-safe observ- 
ables of the (2+l)-jet type. In next-to-leading order, these programs are still restricted 
to the case of photon exchange. A natural extension would be the case of Z exchange, 
and the inclusion of charged- current processes. Moreover, the next-to-leading order 
is not yet available for polarized initial states, and for the pt spectra of identified 
hadrons. 

Concerning the phenomenology of hadronic final states, the main problem is to 
understand the phase space region of small Q 2 (where, for example, the Thrust distri- 
bution is not well described in next-to- leading order), and the region of forward jets. 
With high-statistics data, the extraction of physical quantities can be improved by 
hard cuts to remove these dangerous regions, but a better theoretical understanding 
is certainly desirable. The physics of deeply inelastic hadronic final states will thus 
continue to be a very interesting topic in the future. 
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